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We demonstrate that ultrafast carrier excitation can drastically affect electronic structures and 
induce brief surface plasmonic response in non-plasmonic metals, potentially creating a plasmonic 
switch. Using first-principles molecular dynamics and Kubo-Greenwood formalism for laser-excited 
tungsten we show that carrier heating mobilizes d electrons into collective inter and intraband tran¬ 
sitions leading to a sign flip in the imaginary optical conductivity, activating plasmonic properties 
for the initial non-plasmonic phase. The drive for the optical evolution can be visualized as an 
increasingly damped quasi-resonance at visible frequencies for pumping carriers across a chemical 
potential located in a d-band pseudo-gap with energy-dependent degree of occupation. The subse¬ 
quent evolution of optical indices for the excited material is confirmed by time-resolved ultrafast 
ellipsometry. The large optical tunability extends the existence spectral domain of surface plasmons 
in ranges typically claimed in laser self-organized nanostructuring. Non-equilibrium heating is thus 
a strong factor for engineering optical control of evanescent excitation waves, particularly important 
in laser nanostructuring strategies. 


I. INTRODUCTION 

Collective non-equilibrium effects on transport prop¬ 
erties in metallic systems are topics of current inter¬ 
est as non-standard behaviors play an increasing role in 
optics, thermodynamics or magnetism. Ultrafast laser 
electronic heating can mobilize localized states, induc¬ 
ing massive nonlinear contributions to optical and ther¬ 
mal transport or structural metastability. Understanding 
the intimate electronic mechanisms of optical coupling 
in excited solids and controlling excitation transients is 
therefore of prime importance as new applications are 
emerging in ultrafast optics, plasmonics and heat control 
strategies. If consequences for transport are readily ex¬ 
pectable [Ei, and the role of nonequilibrium electronic 
transfer was early recognized Si , subtle atomistic ef¬ 
fects can be inferred at the level of the band structure 
with unexpectedly strong consequences. Charge trans¬ 
fer and screening during non-equilibrium concur to a dy¬ 
namic self-adjustment of thermal and optical properties 
to accommodate swift excitation. Recoules et al. S 
proved that participation of localized d orbitals in no¬ 
ble metals enforces the lattice cohesion. Non-equilibrium 
charge supply allows large variations of thermal charac¬ 
teristics in transition metals Q , fluidifying energy trans¬ 
port Q. Non-thermal distributions can restrict the colli- 
sional phase-space around the Fermi level, severely damp¬ 
ing electron-phonon coupling Q. Furthermore, electronic 
occupation of delocalized states and filling from localized 
reservoirs (e.g. d-bands in transition metals) redefine 
classical views on free electron density introducing 
significant dynamics in optical behaviors. Tuning ultra¬ 
fast optical response via electronic reactions was thus 
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proposed for active large bandwidth modulation in ul¬ 
trafast plasmonics [T^ . Engineering macroscopic optical 
response on ultrafast scales equally relates to the onset 
of optical resonances on excited surfaces with 

spectral and spatial disturbances, e.g. Wood’s anomalies. 
One debated example concerns ultrafast laser nanostruc¬ 
turing of solids in self-assembled regular patterns. The 
universal phenomenon - observed half a century ago 11311 
- carries application potential in functional surfaces |14| , 
color-coded optical traceability and multidimensional in¬ 
formation storage [ilH, or feedback-driven nanolithog¬ 
raphy [17|. Key questions in laser-induced periodic sur¬ 
face structures (LIPSS) concern the origin of spatially- 
modulated energy patterns [l8| and the interplay be¬ 
tween optical resonances and nonlinear feedback. Elec¬ 
tromagnetic calculations indicate the involvement of sur¬ 
face waves inducing collective motion [T^ . Among them, 
surface plasmon (SP) with its ability to enhance electro¬ 
magnetic energy on nanoscales holds a determinant role. 


The potential involvement of SPs interference in laser 
nanostructuring is in close relation to surface opti¬ 
cal indices. Thus the experimentally-observed onset 
of polarization-dependent periodic structures for non- 
plasmonic metals, as demonstrated on tungsten M , is 
intriguing in optical ranges theoretically forbidden by 
predictions based on ambient optical indices. A typ¬ 
ical example of ultrafast LIPSS at 800 nm is provided 
in Fig. [IJa), showing quasi-wavelength periodicity. The 
positive real permittivity of W precludes normally the 
generation of electronic surface waves and plasmon po- 
lariton coupling. The possibility to excite collective ma¬ 
terial motion pinpoints to a more intricated dynamic re¬ 
sponse of the system. We investigate the impact of ul¬ 
trafast heating on electronic structures capable of ini¬ 
tiating a large excursion of optical properties, namely 
an excitation-driven plasmonic state in otherwise non- 
plasmonic solid non-Drude-like metal. First-principles 
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approaches are chosen to interrogate electronic-driven 
evolution of optical properties and their consequences on 
plasmonic behaviors at time and spatial scales difficultly 
accessible by experiments. 


II. CALCULATION DETAILS 

The calculations correspond to ultrashort pulse ir¬ 
radiation conditions enabling LIPSS around the dam¬ 
age threshold (90 mJ cm“^ absorbed fluence). Tran¬ 
sient electron and lattice temperatures in W for sin¬ 
gle ultrashort pulse irradiation (Fig.[TKb)) are estimated 
using a two-temperature hydrodynamic code (Esther) 
( 2 ^ which describes the energy balance using Helmholtz 
optical formalism and electron-phonon relaxation with 
temperature-dependent transport properties [1, . Elec¬ 

tronic temperatures on the rising heat cycle can amount 
to 2.5 X 10^ K, inducing electron-phonon nonequilibrium 
up to 3 ps, with the material in solid state for at least 1 ps. 
Internal electronic thermalization is assumed, justified 
by high levels of energy deposition in the vicinity of the 
threshold which accelerates electronic energy exchange to 
fs scales In this range, calculations are carried out 
with the plane-wave code Abinit [2l| , in the frame of the 
density functional th eory extended to finite elec¬ 

tronic temperatures |24l |. The generalized gradient ap¬ 
proximation is employed to model exchange and cor¬ 
relation energies, and projector augmented-waves atomic 
data [ 2 ^ account for the effects of nuclei and core elec¬ 
trons. A 54 atoms supercell of body centered cubic (bcc) 
W is considered, with 5s^5p^5d^6s^ valence electronic 
configuration; 4/ electrons are n^lected as their effect is 
weak at the Tg considered here Q. We first perform ab 
initio molecular dynamics simulations in the isokinetic 
ensemble at room temperature, with Ti = Tg = 300 K 
during 2 ps. This calculation provides an average ther¬ 
modynamic equilibrium from which ionic configurations 
are extracted, standing for representative states of the W 
lattice at ambient conditions. Then electronic structures 
are computed at Tg of 300, 10"^ and 2.5 x 10^ K. They 
serve as a basis to evaluate optical properties in solid 
phase at these levels of electronic excitation. 


III. RESULTS 


relates to the roughly symmetric profile of the DOS on 
both sides of the electronic chemical potential p. The 
high density of empty electronic states on the right side 
of the Eermi level collects the excited electrons from the 
left-sided filled electronic states, leading to a weak Tg 
dependence of p, contrary to most of transition metals 
0. The location of the chemical potential within a 
mid-range pseudo-gap between high density of filled and 
empty electronic states determines the non-plasmonic be¬ 
havior in the 1-5 eV range, i.e. for visible photon energies 
bridging the gap. The increase of Tg determines how¬ 
ever a slight augmentation of the number of d electrons 
from low-lying part of sp orbitals, implying a moderately 
strengthen localization of the charge density, noted in in¬ 
creasing Hartree energies @ . The main Tg effect resumes 
to a Fermi broadening within the DOS. The impact of the 
occupation probability is discussed below. 
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FIG. 1. (Color online) (a) Typical field-perpendicnlar snrface 
periodic nanostructures at quasi-wavelength periodicity in¬ 
duced by 800 nm 50 fs laser pulses on W. (b) Electron (dashed 
red) and lattice (solid black) temperature evolutions for W 
irradiated by a 50 fs laser pulse. Dotted line indicates the 
standard vaporization limit, (c) DOS of bcc W computed 
at Ti = 300K and its evolution with Tg. Dashed lines indi¬ 
cate the Fermi levels and dark colored areas show occupied 
electronic states mainly impacted by 800 nm photons. 


We first focus on the evolution of the electronic den¬ 
sity of states (DOS) with electronic temperatures in the 
presence of a cold, ambient lattice, with calculated DOS 
profiles provided in Fig. [T](c). If finite ionic tempera¬ 
ture induces a certain smoothing of the features com¬ 
pared to calculations at T^ = OK @ due to a loss of 
degeneracy arising from atom oscillations around their 
high symmetry position, a highly-structured d band re¬ 
mains. The DOS shape is remarkably stable against elec¬ 
tronic heating with a maximal d block shift of 0.15 eV at 
2.5 X 10“^ K, indicative of a stiff electronic structure. This 


Relying on the as determined electronic structures, the 
subsequent optical properties are obtained from an av¬ 
erage on three ionic configurations. The real part of 
the frequency-dependent optical conductivity is obtained 
within the Kubo-Greenwood (KG) formalism 
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where electronic transitions from i to j states are inte¬ 
grated over the reciprocal space for each photon energy 
huj, accounting for the Fermi-Dirac occupations /(e) and 
the eigenenergies of the electronic states e. V repre¬ 
sents the velocity operator, il is the volume of the cell 
while a corresponds to the three spatial dimensions. The 
imaginary part of the frequency-dependent conductiv¬ 
ity is obtained using the Kramers-Kronig (KK) relation, 

ai{uj) = /q°° dw', where V is the principal 

value of the integral. Frequency-dependent permittivi¬ 
ties and optical indices can be derived. 

Optical properties follow electronic structure evolu¬ 
tions and large excursions up to plasmonic states are 
argued below. Fig. [5] shows the computed frequency- 
dependent optical conductivities a = ar + i(Ji and optical 
indices n = n+ik as a function of the electronic tempera¬ 
ture. A good agreement is found between theoretical val¬ 
ues obtained at 300 K and the reflectivity measurements 
of Weaver [2^ , confirming a realistic description of the W 
electronic structure, ar is an implicit measure of optical 
absorption and its spectral behavior can be inferred based 
on the profile of the energy bands calculated in Fig.[Tl(c). 
For the low Te case (dotted blue line), at photon energies 
below 0.5 eV, the intraband part dominates, resulting 
from electronic transitions inside 6sp bands. This contri¬ 
bution rapidly decays as the photon energy increases. On 
the other hand, interband transitions originating from 
the partially filled 5d sub-bands become gradually more 
important [^. From 0.5 to 5.1 eV, photon-driven elec¬ 
tronic transitions access an increasingly larger DOS do¬ 
main (±/ii^) centered on /r, resulting in a stepwise build¬ 
up of C7r(uj) mapping the local DOS. This trend can be 
seen as the wing slope of a resonant behavior around the 
5.1 eV peak given by the particular electronic DOS split¬ 
ting around the Fermi level and the finite width of the 
sidebands. The conductivity finally decreases once main 
d-bands peaks have been included and diluted into the 
continuously increasing electronic transition phase-space. 
From the KK relation, cri(uj) relies on the overall ar(uj) 
profile. Conceptually, this integral can be separated in 
two parts depending on the value of oj' with respect to 
the reference value w. For co' € [0,a;[, the integral is neg¬ 
ative and reverse signs for higher uj'. Accordingly, the 
magnitude of ai depends on the profile of the real con¬ 
ductivity mainly around uj. A quasi-symmetric profile of 
ar tends to balance positive and negative components of 
the integral and provide low values for ai. On the con¬ 
trary, an asymmetric profile emphasizes the respective 
sign components of the integral, depending on the trend 
of asymmetry (left or right-turned), determining positive 
or negative values of the imaginary conductivity. For W, 
the positive slope of ar profile from 1 to 5.1 eV with val¬ 
ues ranging from 0.5 x 10^ to 1.7 x 10"^ S cm“^ leads to 
a negative ai(uj) at ambient conditions. A spectral view 
on ai(uj) indicates an anomalous-like dispersive behavior 
related to the absorption resonances. 

The profiles of optical conductivities in Fig.[2Ka,b) and 
optical indices in Fig. [lKc,d) significantly change with 
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FIG. 2. (Color online) Real (a) and imaginary (b) part of 
the frequency-dependent optical conductivity; real (c) and 
imaginary (d) part of the optical indices. Colored curves 
stand for the electronic temperatures of 300 K, 10^ K and 
2.5 X 10“^ K and the squares represent the experimental data 
for non-excited W [^ . 


the electronic temperature, requiring an introspection in 
the corresponding excitation-driven electronic effects. If 
the observed changes in optical properties are not a con¬ 
sequence of DOS variation due to its stability against 
heating, the electronic temperature affects the degree of 
filling. Thus the optical evolution is the direct conse¬ 
quence of a broaden Fermi-Dirac charge redistribution, 
making more states available for photon-induced transi¬ 
tions. With the increase of Te, the optical transitions 
concern an increasing /r-centered interval mainly corre¬ 
sponding to [—hv — ^kTe, hv IfcTe] [s^, illustrated 
in Fig. [IKc) for a photon energy of 1.55 eV. The conse¬ 
quence is manyfold, observable in the cr^ which maps the 
spectral absorption. Firstly this produces a broadening 
of available transition phase-space for the various photon 
energies. This translates into an increase of the intraband 
part, balanced by a decrease in the interband domain 
(qualitatively similar to a clockwise turn), damping the 
ar 5.1 eV peak. The subsequent asymmetric change of 
slope in ar with Te determines the increase of ai in the 
1-5 eV (as the dispersive behavior across the resonance 
flattens) and a passage in the positive values domain for 
a significant part of the spectral domain. The situation 
becomes more clear in the evolution of the optical in¬ 
dices. At 1.55 eV (800 nm), relevant for the LIPSS, op¬ 
tical indices are strongly affected by the heating of the 
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electronic system, with n = 4.0 + 2Ai at 3 x 10^ K reach¬ 
ing n = 2.1 -I- 3.8* at 2.5 x 10^ K. Thus, in the near 
infrared and low frequency visible part of the spectrum, 
the real part of the index goes down and the imaginary 
part augments. 

IV. EXPERIMENTAL CONFIRMATION 

The evolution was confirmed by time-resolved pump- 
probe ultrafast ellipsometry on laser-excited W surfaces. 
The transient optical properties upon ultrafast (120 fs) 
laser irradiation were probed at 1.55 eV photon energy 
using a two-angle one color time-resolved ellipsom etry 
method following the technique proposed in Ref. [^. 
The static properties of the non-excited surface were 
first evaluated ex-situ using a commercial ellipsome- 
ter (Uvisel, Horiba Jobin Yvon) and the results give 
h = 3.57 -I- *3.15 for massive W materials (Goodfel- 
low 99.95% purity, mechanically polished). Alongside 
massive W, foils (Goodfellow 0.3 mm thick, polished at 
0.1 mRa -purity 99.95%) were equally used as reference. 
The dynamic reflectivity changes were subsequently in¬ 
terrogated on massive thick targets by p-polarized 120 fs, 
800 nm low energy non-perturbing probe laser pulses at 
27.1° and 65.8° incidence angles. The probe pulses 
were time-synchronized with fs accuracy with the ex¬ 
citing pump pulse of equally 120 fs 800 nm, arriving at 
normal incidence on the W surface. Two photodiode de¬ 
tectors were used in imaging geometries with respect to 
the surface. The probed zone is significantly smaller than 
the spatial extent of the excited region. The exciting flu- 
ences were chosen slightly below the ablation threshold. 
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FIG. 3. (Color online) Example of dynamic time-resolved re¬ 
flectivity traces at the two probe angles (27.1° and 65.8°) 
at incident pump peak fluence of 0.12 J/cm^. Inset: (n, fe) 
contour plots corresponding to measured reflectivities at the 
two given angles. The intersection point represents a uniquely 
determined (n, k) pair satisfying simultaneously the two re¬ 
flectivity conditions. Non-excited values obtained from ellip- 
sometric measurements on bulk and foil samples as well as 
literature data are given as references. 


The swift optical activity on massive W samples is 
underlined by strong reflectivity changes, with a reflec¬ 
tivity snapshot example for an input peak fluence of 
0.12 J/cm^, just below damage threshold being shown in 
the Fig. [3] (a). The time-profile of the reflectivity tran¬ 
sients maps the heating/cooling cycle of the electronic 
system and the associated redistribution in the electronic 
occupation around the chemical potential with T^. The 
measured maximum values of the transient reflectivity 
changes occurring just after the peak of the excitation 
pulse were then depicted in terms of corresponding (n, k) 
values in the (n, k) optical phase-space. They represent 
uniquely-determined (n, k) pairs obtained by inverting 
Fresnel formulas at the given angles. The result is pre¬ 
sented in Fig. 13] (b). The intersection point of the (n, k) 
contour plots corresponding to two angle-resolved mea¬ 
sured reflectivities were used to extract the correspond¬ 
ing (n, k) pair, allowing thus access to the evolution of 
optical indices. The result is h = 3.21 -|- *3.77, where 
resulting optical indices indicate slight n decrease and 
increase in k above the experimental value of its disper¬ 
sive part. Note that these correspond to a macroscopic 
state averaged over a region set by the optical penetration 
depth of 18 nm, with inhomogeneous temperature distri¬ 
bution. The accuracy of the measurement is affected by 
the roughness and local planarity of the surface and care 
was taken to minimize the errors. The optical evolu¬ 
tion has strong impact on the possibility to excite sur¬ 
face plasmon, with the fulfilment of the required optical 
conditions. 



FIG. 4. (Golor online) Real part of the effective refractive in¬ 
dices (rj) in case of air-tungsten interface as a function of Te. 
The symbol line stands for experimental data. Inset: varia¬ 
tion of plasmon wavelength as a function of laser wavelength 
in a spectral domain ranging from 0.9 to 2.3 eV. 


For air-material interfaces the condition for surface 
plasmon existence reduces to Er = — k'^ < —1, a. con¬ 

dition fulfilled in the conditions of Fig. [31 marking thus 
the ultrafast plasmonic activation at fluencies in the close 
vicinity of the damage threshold. 
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V. DISCUSSION 

Surface plasmon periodicity Asp as a function of laser 
wavelength A is given by Asp = A/ry, where rj = 
5R{[n^/(h^ + 1 ) 11 /^} is the real part of the effective re¬ 
fractive index [^. Calculated 77 is plotted in Fig. |4] and 
the absence of data indicates non-existence domains of 
surface plasmons. With the increase of Tg, one can note 
the expansion of the existence domain in the red photon 
energy region. The temperature-induced broadening of 
the electronic transition domain determines thus a plas- 
monic switch at visible optical frequencies. It appears 
that, upon electronic excitation, the light-induced onset 
of plasmonic character can sustain an origin based on op¬ 
tical resonances for LIPSS at 800 nm even though room 
temperature optical indices indicate non-plasmonic prop¬ 
erties. As expected for air-metal interfaces, 77 remains 
close to one, especially at low photon energy, leading 
to Asp slightly inferior to the laser wavelength (inset of 

Fig. 15 ■ 

From the expression of the dielectric permittivity e = 
l+ia/u!£o the existence condition can be expressed more 
directly in terms of imaginary conductivity (and implic¬ 
itly on real permittivity) with ai > 2 ui£q. To clarify pro¬ 
cesses leading to the increase of the existence domain, we 
propose a simplifying approach. The imaginary conduc¬ 
tivity is tentatively split into intra and interband com¬ 
ponents as depicted in Fig. [5Ka,b). We first compute 
optical properties for a degenerate electronic system at 
T, = OK and = 300, 10^ and 2.5 x lO^K. This dis¬ 
regards ionic temperature and phonon effects, assuming 
they assist mostly the intraband component via momen¬ 
tum conservation conditions. Conceptually this is justi¬ 
fied in a classical view by the temperature dependence of 
damping frequencies. Accordingly, the intraband contri¬ 
bution vanishes and ~ cr,(r, = OK) (Fig. [SKb)). 

The intraband part at Ti = 300 K is then extracted 
by subtracting the as-determined T^-insensitive interband 
part from the total value of the imaginary conductivity 
(Fig. nib)). In the visible spectral range, the intraband 
part rapidly increases with the rise of Tg and sat¬ 

urates at 10^ K. An example of this behavior at 1.55 eV 
is given in the inset. At the opposite, the interband part 
shows a more complex behavior, with negative values of 
^mter 0 to 5 eV energy interval that leads to nega¬ 

tive values of the total imaginary conductivity, explaining 
the non-plasmonic nature of this metal. At higher en¬ 
ergies, a resonance-like dispersive shape induced by the 
5.1 eV peak of is clearly visible (Fig. [Sib)). By weak¬ 
ening the peak of Ur, Fermi broadening also flattens the 
shape of cr?"*®’’ around the resonance and gradually re¬ 
duces its negative component, progressively switching on 
plasmonic properties of W. The increase is not yet satu¬ 
rated at 2.5 X lO'^K giving the possibility of a stronger 
lift and a subsequent extension of the plasmon existence 
domain for higher electronic temperatures. This high¬ 
lights the preponderant role of d-electron-driven inter¬ 
band transitions in the optical response at high T^. 
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FIG. 5. (Color online) Intra and interband contribution to 
the imaginary conductivity as a function of K. Inset: Local 
evolution of ai{Te) at 1.55eV. 


The electronic evolution described above is related to 
the fast achievement of the non-equilibrium electron- 
lattice phase. However, W is a strong coupling material 
and lattice effects may be rapidly generating along the 
electronic relaxation. Performing simulations at different 
lattice temperature values we observe that a ionic tem¬ 
perature effect is indeed observable, increasing slightly 
the plasmon existence domain. This is attributed to a 
loosening of the atomic order initially responsible for the 
material non-plasmonic nature. However the effect is of 
secondary importance compared to the charge redistri¬ 
bution induced by the electronic temperature within the 
band structure. We conclude that lattice heating, al¬ 
beit its influence on atomic order, does not severely im¬ 
pact the calculated electronic structure on the relevant 
timescales, maintaining a dominant electronic drive for 
the process. 


VI. CONCLUSION 

In conclusion, electron temperature dependent ab ini¬ 
tio MD-KG calculations of solid state W indicate large 
duty-cycle optical tuning upon ultrafast laser irradia¬ 
tion and potential excitation of propagating collective 
electronic motion for an initially non-plasmonic state. 
We demonstrate that the evolution of optical properties 
with Te produces, via the necessary optical conditions, 
an extension of the predictable existence domain of sur¬ 
face plasmon in the visible range, rendering possible a 
transient plasmonic phase and a potential plasmonic in¬ 
volvement in LIPSS. The necessary condition for opti¬ 
cal indices was confirmed by time-resolved ellipsometry 
on excited W surfaces. The dynamic evolution mecha¬ 
nism is related to a redistribution of localized d-electrons 
across a chemical potential located in a d-band pseudo¬ 
gap. Since the DOS is not distorted by electronic heat¬ 
ing, the changes of plasmon properties, indicated here by 
corresponding changes of optical conductivities, are prin- 
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cipally due to Fermi broadening within a structured d- 
block, extending the transition space for visible frequen¬ 
cies. By depopulating low-lying electronic states in favor 
of high electronic states, quasi-resonant transitions from 
regions of high DOS are diluted in a continuum of tran¬ 
sitions. The dissimilar behaviors of intra and interband 
absorption events mark thus the evolution of the optical 
conductivity towards fulfilling resonant conditions. Simi¬ 
lar phenomena may occur in other metals exhibiting non- 
plasmonic characteristics [s^, especially Cr and Mo, as 
they crystallize in a similar structure. Stronger changes 
of optical properties [sj are to be expected in case of 
transition and noble metals where the electronic struc¬ 
ture varies substantially with electronic temperature. All 
these effects have implications for optical tunability and 


ultrafast switching of plasinonic properties extending be¬ 
yond the domain of laser nanostructuring and validate 
concepts of structure engineering for metallic materials. 
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